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Abstract 

Starting from a phenomenological Hamiltonian originally written 
in terms of angular momentum operators we derive a new quantum 
angle-based Hamiltonian that allows for a discussion on the quantum 
spin tunneling. The study of the applicability of the present approach, 
carried out in calculations with a soluble quasi-spin model, shows that 
we are allowed to use our method in the description of physical sys- 
tems such as the Mnl2-acetate molecule, as well as the octanuclear 
iron cluster, Fe8, in a reliable way. With the present description the 
interpretation of the spin tunneling is seen to be direct, the spectra 
and energy barriers of those systems are obtained, and it is shown 
that they agree with the experimental ones. 
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1 Introduction 



Besides the well-known examples of tunneling processes that permeate the 
literature in molecular and nuclear physics PQ 12 E] , in recent years it has been 
shown that quantum tunneling also appears in as broad physical situations 
as magnetic moments in spin glass systems, single magnetic impurities in a 
crystal field, a single domain, a ferromagnetic, or antiferromagnetic, domain 
wall, and also in spin tunneling in molecules IH El E] , emphasizing the im- 
portance of the angle/angular momentum degree of freedom in these cases. 
This quantum effect occuring in such mesoscopic, and even in macroscopic 
systems points to the necessity of a discussion on how such tunneling can be 
realized and controlled 0. 

In what concerns the case of spin tunneling, or the quantum tunneling 
of magnetization, QTM, the theoretical approaches have pointed to some 
different ways of treating this problem. In particular, the spin tunneling was 
first treated by van Hemmen and Siito jH], and by Enz and Schilling 
From different perspectives, various approaches were introduced: by using a 
generalization of the usual WKB method adapted to spin systems llUllTH IT^. 
or by using Feyman's path integral treatment of quantum mechanics and 
also by using coherent states Jl]. On the other hand, from the experimental 
point of view, great advances have been made since the synthesis of the 
Mnl2-acetate molecule in 1980 ^Hl- Among the properties of this molecule 
it emerged an indication of spin tunneling that has strongly motivated great 
efforts in characterizing and elucidating the pure quantum contribution to 
that phenomenon ITBl El ^1 El I2D1 • The discovery of the same indications of 
spin tunneling in the octanuclear iron cluster, Fe8 [2I1I221, has thus reinforced 
the proposed idea of pure quantum spin tunneling, and formal treatments 
involving phenomenological Hamiltonians were then widely proposed and 
discussed 

In any case, the angular momentum operators obeying the standard com- 
mutation relations 

constitute the fundamental blocks from which the starting Hamiltonian de- 
scribing the spin system is constructed. In general, a Hamiltonian is proposed 
based on experimental considerations on the symmetries of the physical sys- 
tem, and also under the assumption of low temperature such that the pure 
quantum effects dominate. In this sense, a general Hamiltonian encompass- 
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ing the molecule symmetries plus possible paralel and tranverse magnetic 
fields may be written as 

H = Ho + CiJ, + C2Jx, (1) 

being Hq the molecule Hamiltonian, and the constants Ci measure the in- 
tensity of the additional external magnetic fields. This description clearly 
separates and takes into account the molecule total angular momentum de- 
gree of freedom as the dominant one and, as far as no other degree of freedom 
is directly and explicitly involved, it is considered a sound description of the 
system. As such, all of its kinematical content can be described in terms of 
combinations of the set of 2j + 1 states generated by all the projections of the 
total angular momentum operator J, which obeys the eigenvalue equation 

j2|jm) = fi^j(^j + l)\jm) . 

It is then clear that those Hamiltonians can then be diagonalized within the 
state space defined by { \jm) } whose dimension is given by the total angular 
momentum J. At the same time, as it is well known, since the Hamiltonian 
is written in terms of the generators of that algebra, it commutes with the 
Casimir operator thus giving at least one more constant of the motion - 
besides the energy; in fact, all the symmetries present in the Hamiltonian 
must also play a precise role, giving rise to characteristic features in the 
energy spectrum. 

In this paper we want to extend the content of previous publications [211 
that proposed to construct and present the main features of a potential 
function, described in terms of an angle variable, which can be used to picture 
in a more intuitive way the spin tunneling. In those works, we have shown 
how to obtain such a potential and how pure quantum tunneling shows up 
as a transmition through an energy barrier. In the present contribution we 
show how to extract a new Hamiltonian which also includes a kinetic term 
- not only a potential energy function - also described in terms of an angle 
variable, from a phenomenological Hamiltonian, Eq.(PQ), that can be used to 
study the spin tunneling. 

To construct such a Hamiltonian, as for the kinematics, we use the su{2) 
coherent states since they constitute an overcomplete basis of states, 
and also because they present the desired properties for the study of the sys- 
tems of interest in the semiclassical limit. At the same time they also include 
quantum effects such as zero point energies associated with the corresponding 
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angle and angular momentum variables. We also use the Generator Coor- 
dinate Method (GCM) of wide use in nuclear and molecular physics |?71 OH] 
mainly because it has been shown to be a useful tool when one deals with 
collective behavior in many-body systems, that is the case in the spin systems 
of interest where the tunneling degree of freedom has a suggested collective 
character. 

From the operational point of view, the GCM is used so that an energy 
surface is generated that embodies off-diagonal matrix elements that account 
for the quantum corrections to the diagonal contribution which, by its turn, is 
basically dominated by the semiclassical contribution. The present approach 
proposes to take into account those off-diagonal contributions by calculat- 
ing the moments of that energy distribution; the zeroth-order moment is 
associated with the potential function, while the second-order moment is as- 
sociated with an inertia function: together they constitute our approximate 
Hamiltonian. It is clear that the higher order terms can also be obtained 
so that the complete series gives rise to a new representation of the starting 
phenomenological Hamiltonian. A truncated series up to the second term is 
then an approximate description of the system, and its validity must neces- 
sarily be discussed. It is precisely the use of the su{2) coherent states that 
makes this study feasible since our moment expansion of the energy surface 
is then akin to the expansion discussed in the literature which also 
makes explicit use of the coherent states. In the present context, it is the 
total spin of the molecule, j, the quantum number that determines how many 
terms we must keep in the moment expansion series in order to get a reliable 
Hamiltonian. 

It is worth noticing that the present approach is based on quantum 
grounds from the beginning in such a way that it does not need to be 
quantized in any form. Therefore, it has a different perspective from the 
approaches that start from a classical Hamiltonian and then quantize it. 
Furthermore, what is important to keep in mind in the present context is 
that the associated pair of quantum operators related to the rotation degree 
of freedom must be treated properly [201 EH E2]- 

This paper is organized as follows. In Section II we introduce the su{2) 
coherent states in the discussion of the semiclassical approach for a general 
phenomenological Hamiltonian, and show how some preliminary results are 
obtained which suggest that quantum corrections may be important in some 
cases. The GCM is presented in Section HI where we show how the energy 
and overlap kernels of the method are obtained for the general phenomenolog- 
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ical Hamiltonian, and the new general angle-based Hamiltonian is obtained. 
In this section we also carry out the calculations for the soluble Lipkin quasi- 
spin model [331 El] since it is a valuable testing ground for systems of the kind 
we are interested in, and also for establishing a reliable truncation criterion 
for the moment series in terms of the value of the spin. The result that we 
obtain guarantees the application, in section IV, of the second-order approx- 
imate Hamiltonian to the description of two interesting molecules studied in 
the literature, namely, the Mnl2-acetate, and the Fe8 cluster, and we show 
that we reproduce the experimental results to a good degree of accuracy. 
The spin tunneling phenomenon occuring in these molecules is qualitatively 
discussed in this context, and the main aspects of the potential barrier are 
discussed. The role of the effective inertia function in the tunnehng process 
is also stressed. Finally the conclusions are presented in section V. 

2 A semiclassical approach 

As a starting point to get a semiclassical approach for describing spin systems 
we may consider the su{2) coherent states defined in a general form as j2ni 

\jz) = {l + \zfy'e-'^\j, -j), 

or also as 

where z is a complex variable given by 

z = tan — e~*^, 
2 

and j characterizes the angular momentum state multiplet. These states are 
nonorthogonal and normalized, and furthermore satisfy an overcompleteness 
relation 

j \jz){jz\dfx{z) = 1, 

which allows us to use them as an useful state basis. 

Moreover, these states are a very interesting starting point when semiclas- 
sical aspects of the physical system must be discussed, as has been already 
discussed in the literature PU] . 
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Let us consider now a spin Hamiltonian of the type that, for simphcity, 
does not include transversal magnetic fields 



H=AJ, + BJl + G{Jl + J!) (2) 
written in terms of the basic operators obeying the su{2) algebra 

[J„ J±] =±J± 

and 

[J+, J_] = 2J,. 

It is clearly seen that by properly choosing the parameters of that Hamilto- 
nian we end up with particular expressions associated with physical systems 
we may want to study, for instance if we choose i? = we have the Lipkin 
quasi-spin solvable model 



Hl = AJ, + G[Jl + Jl) . 

By choosing G = we get the Mn-12 acetate model, where a paralel magnetic 
fields is present jj], namely 

Hac = + bjI , 

and we also see that by choosing A = we have a model Hamiltonian for 
the octanuclear iron cluster, [Fe802 (0H)j^2 (tacn)g]®+, or just Fe8 |2Ill22]j 
in the absence of external magnetic fields 

Hpe = BJl + G (4 + Jl) . 

We then recall that a description based on an angle representation of these 
models, out of that general Hamiltonian, can be obtained by directly calcu- 
lating its respective normalized matrix elements with the proposed coherent 
states. To this end, we collect the basic results 

{3AJz\3z) 
/ . I . \ = -J cosa , 

{jz\Jz\jz) .2 ■ f ■ IV 2 



{jz\jz) V 2 
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{jz\Jl\jz) 



a 



a 



and 



{3z\3z) 

{3z\Jl\]z) 
{jz\jz) 



2j (2j - 1) sin - cos -e 



{jz\Jl\jz) 
{jz\jz) 



Using these results we get the energy surface associated with the Hamilto- 
nian, namely, 



—Aj cos a + B 
2G2j (2j - 1) sin 



f-J 



J 



sin^ a 



+ 



(3) 



2 a 2 



— cos — COS 2f . 
2 2^ 



At this point, we see that we can determine the value of ^ that minimizes 
the energy surface, leading to an energy curve in the variable a. 

Let us analyse this expression in a particular case. It is evident that, if 
we choose the parameters so that we describe the Mnl2-acetate molecule in 
the absence of paralel and transverse magnetic fields, i.e., A = G = (the 
energy surface in this case is ^-independent), the energy curve 



n(a) 



B 



J -J 



sin^ a 



presents minima at a = 0, vr for 5 < 0. For the values B/ks = —0.6 K, 
where is the Boltzmann constant, and j = S* = 10, as found in the 
literature ^^1; we obtain the energy minima at -Emin = —60.0 K. This is, 
however, precisely the ground state energy value obtained by diagonalizing 
the Hamiltonian associated with the corresponding model, Eq.(0), within the 
|j(= S)m) states basis. This result thus indicates that such a semiclassical 
approach does not embody quantum content enough to describe such system 
because, as we see, the bottom of the energy curve coincides with the ground 
state energy. This fact shows that the quantum correlations related to the 
rotation degree of freedom may be important and, furthermore, it is not 
taken into account in this approach. 



3 Formalism 

Let us now go one step further in order to construct an energy surface that de- 
scribes the spin system and at the same time that encompasses more quantum 
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correlations than the previous procedure does. In this sense the Generator 
Coordinate Method (GCM) [27] extends that approach because it also em- 
bodies off-diagonal elements, i.e., if we start from the su{2) coherent states 
we have to calculate the so-called energy and overlap kernels, as already 
discussed [ 



N{a,a') = {jz'\jz) = {ja'^lja^ = 

respectively. It must be observed that hereafter we will only consider a as our 
generator coordinate. The variable ^ will be properly chosen at the beginning 
in each case so as to minimize the energy surface as pointed out before. 

In the same way as we did in the last section, the basic matrix elements 
of the operators are collected 

/ ■ /I T I ■ \ ■ 2i-i fot' - a\ f a' + a\ 

\JZ \ Jz\]z) = -J COS I 1 COS I — - — 1 , 

{jz\JAjz) = - cos ^ j+-{2j -1) cos ^ 



and 



{jzVl + Jl\jz) = 2j (2j - 1) cos'^-^ 



a' ~ a 



2 



2ii ■ 2 2 , -2i£ -2^ 2^ 

e ^ sm — COS — he ^ sm — cos — 
2 2 2 2 

Therefore, the GCM energy and overlap kernels for the general spin Hamil- 
tonian are respectively 

K(a,a'-X) = -Ajcos^'- cos (i^^) + b| cos^. (^^l^) + 

' • - ■;^^J ^ 1 1 O/^G 



B'- (2j - 1) cos^^-^ ( ^ ) cos^ ( ^ ) + 



2Gj (2j - 1) cos^^-^ ( sin' j cos' | + e'^^^ sin^ | cos^ ^ 



(4) 
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and 

N (a, a') = cos^^ f ) • (5) 



The contribution of the off-diagonal matrix elements is immediately recog- 
nized in these expressions, and it can be easily verified that they are less 
and less important as j increases, and that the semiclassical results are thus 
reobtained for a' = a, or when j ^ oo. Thus, the off-diagonal contributions 
may play an essential role in the description of spin systems with small and 
medium values of j since they carry additional quantum information when 
compared with the semiclassical approach. 

Instead of solving the GCM integral equation with these kernels, as usu- 
ally is donej2Zl! we will show in what follows how a spin Hamiltonian - written 
in terms of an angle variable - can be extracted from the GCM kernels. 

Once we have the GCM overlap and energy kernels, we will use a new 
represention for describing the spin system^35]. Following the scheme already 
presented before [211 123, which makes use of Fourier transforms, we get the 
new matrix representation for the spin system 

rr _ _^ S (j, n, n') 

2 r (^i^ + n+n' \ Y pi+1 _ n+n' \ \"n',n+l + On',n-l) + 
\2 2/V2 2/ 

S (j, n, n') ^ 

2 r (j + ^ + 1) r (j - ^ + 1) 

B S{j,n,n') 

Tr (j I "+"') r (j "+"7 ^ '^""''"+2 + '^n',n-2)\ + 

G S (j, n, n') 

I "+"') r (j "+"7 ~ (^"''"+2 + ^n',n-2)] + 

Gj (2j - 1) ^ , (6) 

where 

S (j, n, n') = (j + n + 1) r (j - n + 1) r (j + n' + l)T (j - n' + 1). 

The diagonalization of this Hamiltonian matrix clearly gives the same re- 
sults as those obtained by directly diagonalizing the matrix generated by 
the starting Hamiltonian, Eq.Q, using the angular momentum \jm) states 
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basis. This result is expected because this new representation is obtained 
by the use of unitary transformations acting on the original one. This also 
assures us that this new representation totally preserves the content of the 
original representation and that it is not an approximate version of the energy 
surface. 

This new representation has also the virtue of allowing a further manip- 
ulation that is an essential step in order to write a Hamiltonian in terms 
of an angle variable for the spin system. Namely, to this end we perform 
new Fourier transformations, followed by a discrete version of the continuous 
Weyl-Wigner transformation jSEI- Again, using the results of Refs. [211123, 
we can get an expression for the energy surface - from the GCM energy ker- 
nel - in terms of a new pair of angle variables, namely (3 and (3' . In fact, we 
introduce the new variables ip = 6 = a' — a, u = (3 — (3' and = 

and also I = n + n' and k = n — n' in order to transform the GCM energy 
kernel into the new energy surface. The summations over / that appear are 
restricted so that they run over even/odd values if k is even/odd. 

Since we already have the GCM energy kernel, Eq.(^, we will consider 
hereafter the particular case of no transverse magnetic field, C = 0, for 
simplicity, and therefore we rewrite it as 



Bl + \G\ji2j-l) 



cos2M?)-2|G|j(2j-l)cos^^-^ f^- 



Aj cos^' 1 (^0 cosip + J (2j + \G\^ cos^ ip, 

where the choices for the variable ^ were already carried out. After a direct 
but tedious calculation we obtain the general expression. 



l=-2j l=^2j ^ 
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27r ' cos20^(^j + 2 + 2 / ' 2 / 2 ' 2/ 2 

(7) 

So far no approximation has been made and the above expression contains 
exactly the same physical content as the original phenomenological Hamilto- 
nian. 



3.1 Discussion 

Before carrying out the calculations in order to obtain the moments of the 
above shown energy surface, we must recall some aspects of the Weyl-Wigner 
transformation techniques developed for continuous degrees of freedom jHT]. 
It can be readily recognized that the power series expansion in the momentum 
p that appears in the defining expression of the Weyl-Wigner transform 

H (g, p) = j e'P("-"')/7 (^^^, x-x'^d{x- x') , 

where = q, and H {q, p) is the quantum phase space representative of 
the Hamiltonian operator, can be rewritten in such a form that the new coef- 
ficients are the moments of the energy surface with respect to the difference 
of the original variables, namely 



M„(^^^^ = J H (^^^^,x-x'^ {x-x'Td{x-x'), n = 0,2,4,... 

is the n-th moment of the energy surface in the variable x — x' . Obviously 
we expect that only the even moments do not vanish because otherwise the 
system would violate time reversal symmetry, or, in other words, terms would 
appear in the Hamiltonian that would be proportional to odd powers of the 
momentum. In fact, it is the even parity of the GCM kernels that warrants 
us a time reversal invariant Hamiltonian. 

In the continuous case, by using this procedure one ends up with a series 
the first two terms of which are dependent on the coordinate variable only, 
and proportional to the squared momentum respectively. In some cases, 
the proportionality coefficient in the second term may result in a function 
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of the coordinate variable or, in other words, an effective mass may come 
out. Therefore, these first two terms can be seen as the potential and kinetic 
energy terms of an associated effective Hamiltonian. When these two terms 
dominate the series expansion, i.e., higher order moments can be neglected - 
they may even be completely absent in certain cases - then the approximate 
two-terms effective Hamiltonian can be useful to describe the physical system 
we started with. 

Now, the same can be seen to be valid in the present case of spin degrees 
of freedom if we consider the proper angular function in respect to which 
the moments of the energy surface must be taken. Since the angle is not the 
direct variable to be used here since it does not obey a standard commutation 
relation with the angular momentum operator [301 1^ 122] ^ in contrast with 
the canonical case of ordinary coordinate and momentum discussed above - 
the moments of the energy surface must be calculated with a periodic function 
of the angle. Thus, we perform the calculation with respect to sin(M) ^T] . 
where u = (3 — (3', as already defined. The zeroth order moment will give 
us the potential function in terms of the angle variable, as has been already 
presented in the past [211 123- Since the zeroth order is given by 

2j 

where Am = we end up with 

^^^^ ^ j[B + 2\G\i2j-l)] ^ B-Q\G\ ^, 



n,k=~j Y 

^^p^j{j + l)cos2<l^, (8) 

where = 2j ' + 1 is the angular momentum state space dimension. 

By its turn, the second moment will give us the expression associated 
with the inertia function. In other words, 

27+1'^ 

7(0) = ^ H{u,(f))sm^u Au 

2j 
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leads to 



= (5 - 6 |G|) + 



A cos 



4j(j + l) 



n,fc=— J 

B + 2\G\ 
4 



n (n + 11 



J (j + 1) 



i {i + i)cos; 



6 



3 {j + 1) 



3 (j + 1) 



It is then a direct matter to see that 

/(0) = 



where M((/)) is the corresponding effective "mass" function associated with 
the spin system. 

We now propose that these first two moments must suffice to write a 
Hamiltonian that can, to an acceptable degree of accuracy, describe systems 
consisting of a suitably great number of spins. In this form, in this approxi- 
mation the Hamiltonian is written as 



d 



+ ^(0) 



d(j) \2M{(f)) d(f)^ 

and we straightforwardly get the eigenvalues as well as the eigenfunctions 
ipk (0) by solving the Schrodinger equation 

H{<p)^k {<P) = Ek4Jk (0) 

by a Fourier analysis. As such, the spin number for which the validity of this 
proposal is assured can be estimated by introducing a model Hamiltonian 
and comparing the results obtained from this Schrodinger equation with the 
exact ones, obtained from a direct diagonalization of that model Hamiltonian 
in a \jm) basis, for a great range of values of the spin quantum number j. 
In what follows we will discuss the validity of that approximation using the 
Lipkin mo del [S^. 

As already pointed out before, the Lipkin quasi-spin model is singled out 
by assuming the parameters: A = e, B = 0, and G = ^ = ~ 2(jv^_i) > where 
Ng = 2j, so that the GCM kernels read 



jcos ^ 

2 ^ 



I X AT, 

COS + COS 
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and 



(l + sin^ 9?) — cos^ 



respectively, if we measure energy in unities of e. 

We obtain the potential function directly from Eq.®, which gives 



+ . 2 

— sm I 



2j 



E 

n,k=—j 



Vj (j + 1) 

n{n + 1) 



COS0X 



(9) 



In the same form we get the expression related to the inertia function 

3x , VTOTT) 



+ 



n,k=—j 



2j - 1 ■ 4A^ 



■ COS (/)X 



iv 



n (n + 1) 



J (j + 1) 



xi (j + 1) 

4(2j-l) 



cos 20 



6 



J(j + 1) 



J (j + 1) 



1 



(10) 



With these expressions we can discuss now the approximations proposed 
before. First, if the number of quasi-spins is such that Ng 3> 1, the summa- 
tions over k can be substituted by integrals and we obtain the expressions 



Vl (0) = sm 0- 



and 

h (0) = 

1 ^ 

— v/j(j + i)cos0 J] (-ly 

n=-j 



n=-] 



n{n + 1) 

J U + 1) ' 



N. 



2^ (l + sin2 0) + 



2n + 5 2n - 3 2n + 1 
14 



77. (n + 1) 



J (j + 1) 



respectively. Furthermore, the remaining summations can be seen to have 
well defined limits for A^^ ^ 1, so that we get 

^.W--^oos^-^(^sm=^. (11) 

and 

2 2y 
^ ~ AW ~ + ''''' ^^^^ 

respectively. These expressions are closely related to those obtained from an 
Adiabatic Time Dependent Hartree-Fock (ATDHF) approach to this model 
[SHI- f^ct; expression that gives the inertia function, Eq. (|T^ . almost 
coincides with the one obtained from that approach (A^^ instead of our A^^ — 1) 
while the expression associated with the potential, namely 



shows that the ATDHF mean field approximation does not include the term 
— |cos(/) — ^ sin^ that appears in our approximation. This difference is 
due to the initial GCM ansatz used in our approach which deals with linear 
superpositions of Slater determinants as the starting quasi-spin wavefunc- 
tions instead of pure determinants as it is the case with ATDHF. A similar 
discussion using the GCM has been also presented in Ref. [SH] . 

Figure 1 depicts the potential functions obtained from expressions (jH)) 
and (jllj) , while Figure 2 shows the functions given by and ()12j) , with the 
changed signal, for x = 1-5 and A^^ = 2, 6, 10 and 20, respectively. From these 
figures it can be verified that the curves tend to superpose as A^^ increases 
being that they are almost indistinguishable already for A^^ = 20, i.e., j = 10. 

With these results we can now check the validity of our approach and 
its limiting version when A^^ ^ 1. To compare the two approaches we first 
obtain the exact results by diagonalizing the energy matrix Hny, Eq.l^, 
and, on the other hand, by solving the Schrodinger equation 



1 d f 1 d , 

+ V{(f)) 



2d(f) \M{<i))d<p 



^(0)=E^(0) (13) 



using the expressions for F(0), Eq. (fTT|) . and M(0), obtained from Eq. lfT^ . 
respectively. From the outset it is clear that the second order moment trunca- 
tion proposed here will not suffice to describe the quasi-spin system for small 
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Ng since in these cases higher order terms are also important. Thus the re- 
quirement of large Ng imposed so that a reliable two-terms Hamiltonian can 
be obtained also results in the approximate expressions, Eqs. flTT|) and (fT^ 
respectively. In this form, from a direct inspection of Figures 1 and 2, we 
know beforehand that our approach must be acceptable for A^^ = 2j > 10. 
This guides us to study the spectra of systems with those quasi-spin number 
directly using expressions (jll|) and (jl2j) as our starting point. Therefore, the 
solutions thus obtained give an way of testing the quality of our approximate 
expressions for the potential and inertia functions. 

A comparison between the approximate and exact energies can be seen in 
Figures 3 and 4. A check of our description is shown in Figure 3 where the 
exact ground state eigenvalues of the quasi-spin systems, with A''^ = 10, 20 
and 40, are compared with those obtained from Eq. lfT^ . using Eqs. (fTT|) and 
fll2|) . as a function of the strength parameter x- Iii Figure 4 we show the 
energy errors in the ground state eigenvalues as a function of the quasi-spin 
number A^^ for a fixed value of the strength parameter, namely, x = l-O- The 
results are almost the same for other values of x- It is clear from Figure 4 
that an error as low as 1% is achieved already for a system with A^^ = 8 
{j = 4) thus confirming our previous assumption, and therefore establishing 
a suitable working limit for the number of quasi-spins of the system when 
we use the present approach. This also means that the truncation of the 
series of moments of the energy surface up to the second term can, already 
for 2j = Ng > 8, account for the dominant part of the dynamical content of 
the quasi-spin system. 

Two interesting features of the Lipkin model can be directly discussed in 
this approach. First, it is immediate to see, observing the expression describ- 
ing the potential function, Eq. ljlip . that the strength parameter, namely 
assumes its critical value Xc ^ the minimum at = turns into a maximum - 
when X = Xc = {Ng + 1) / {Ng + 3). This shows that only when A^^ — oo we 
get Xc — ^ 1, as expected Clearly for values x > Xc {Ng < oo) a barrier 
appears at the origin and tunneling can occur under certain conditions ^3] . 
Second, the possibility of tunneling through the barrier, when it is the case, 
can be blocked if the effective mass presents a divergent behaviour at or near 
the classical turning points, i.e., if II (0) exhibits zeroes there. In this form, 
not only the potential barrier governs the tunneling process, but the inertia 
function also plays an essential role in this matter. Due to the behaviour 
of the inertia some regions in the angle domain may not be acessible; this 
can also be seen as a localization of the spin orientation when this singular 
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behaviour occurs: when it turns out to be impossible for the spin to change 
its orientation - when the effective mass diverges - it cannot tunnel. Consid- 
ering II {(p) for the Lipkin model we see that its zeroes occur at 0o = if for 
X = 0, and that they are pushed to the borders of the angle interval as x 
increases to 1, i.e., 0o = i^r for x = 1- For X > 1 there will be no zeroes at 
all. This result shows us that whenever a barrier is present in the potential 
function - and this occurs for X > Xc ^ there will always be a nonvanishing 
tunneling probability. As a consequence, there will always be tunneling in 
the Lipkin model when a barrier is present. 



4 Applications 

4.1 Mn-12 Acetate molecule 

If we choose now the parameters of the phenomenological Hamiltonian such 
that A = gfiBH\\, B = —D [D > 0), and G = we have the already 
proposed Mnl2-acetate molecule Hamiltonian model j31l2Sl- Starting from 
that Hamiltonian, and using the general expression, Eq.((Zj), we end up with 
the corresponding energy surface 



l=-j l=-j 



D 

— COS 

47r 



We can now obtain the expressions for the potential and inertia functions, 

Vac{<I>) = -Dj{j + 1) cos' (P- 



^cos</>V7(JTI) y e'-U^^'^)Jl-!±±}l 
2j + l ^^^^^ V + 
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and 



M(0) 



-D+ ,y^ll^ cos0Vj(j + l) J2 



4(2j + l) 



AT X 



n,k=—j 



n{n + 1) 



J U + 1) 



D 



cos20j(j + l^ 



6 



J (j + 1: 



1 



J u + 1) 



1 



respectively. Since in this case j = S = 10, and all the summations can be 
directly calculated, we are allowed by the estimates of the previous section 
to directly use here the already discussed approximations such that these 
expressions reduce respectively to 



Vac (0) = —S{S + 1)D cos^ — ■ygfisHw cos ( 



(14) 



and 



Iac (0) = —D — ^qhbHw cos — [3D cos 20, (15) 

where 7 = ^S{S + l) (to within 1.4%), ^ = l/S (to within 0.1%), and 
13 = 1.0 (to within 0.47%). Figures 5 and 6 depict Vac{4>) and /ac (0) for 
some values of H\\ where we have adopted the angle domain (— f , ^) for 
the sake of clarity. In Figure 5 we see that the potential for the Mnl2- 
acetate molecule presents two minima whose energy depths depend on the 
applied field. At the same time the expressions for 1^^ (0) shown in Figure 
6 present a similar behaviour. As it is shown, as the paralel field increases, 
one of the minima gets shallower while the other deepens as expected from 
phenomenological considerations 

If we solve the Schrodinger equation with the obtained expressions for the 
potential and inertia function, we obtain Egg = —60.278 K for the ground 
state energy when H\\ = with D/ks = 0.6 K. This result shows that the 
calculated ground state energy is only 0.46% below the exact energy result 
(—60.0 K), obtained by diagonalizing the corresponding matrix. When 
H\\ 7^ the errors are of the same order of magnitude as that. These results 
confirm our previous considerations, and thus allow us to further discuss the 
Mnl2-acetate molecule properties using our approach. However, instead of 
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extracting results from numerical calculations based on the new Hamilto- 
nian, we will discuss some basic properties of this molecule by analysing the 
potential and inertia function analytic expressions. 

Observing that we can consider ^ = 1/5, and that we can introduce 
the field at which the magnetization of the Mnl2-acetate molecule attains 
saturation i.e., 

_2SD 

J^a — , 

we can further rewrite expressions ()14|) and ()15|) as 

on 

Vac (0) = -S{S + l)D cos^ (p - —S>yS{S + l)Hi\ cos0 (16) 

-"a 

2D 

/4^(0) = -2Dcos20- — i/||cos0. (17) 

Ha 

It can be immediately seen that while Vac (0) oc S"^, Iac (0) oc 5*°, thus show- 
ing the relative dominance of the potential energy term in the approximate 
Hamiltonian. 

Now, let us analyse the expressions for the potential and inertia functions. 
First of all let us extract the essential features of the potential function. It is 
direct to see that since the treatment is of quantum nature, it embodies the 
inherent quantum correlations, and as such, even in the absence of paralel 
magnetic fields, the potential minima are —66.0 K, instead of the semiclassi- 
cal value —60.0 K obtained with D/kB = 0.6 K. In fact, in our calculation, 
it is the ground state energy that occurs at —60.278 meaning that this is 
the energy height of the ground state barrier in this description, and further- 
more that the ground state energy does not coincide with the bottom of the 
potential, due to the uncertainty principle, as it must be. This also points 
to the fact that the quantum correction related to the uncertainty principle 
amounts to about 10% of the potential barrier height. In the particular case 
of no paralel field, the minima occur at (pmin = 0, vr, while the maxima occur 

r,f A — ZL 37r _7r\ 
ah v^max 2 ' 2 \ 2 / ' 

The introduction of a paralel magnetic field, ify , only shifts the maxima, 
the minima being kept unaltered, viz.. 



1 H 



arccos 



^S{S + l)Ha 
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In what refers to the inertia function let us study first of all its expression 
by looking for the extrema. Let us consider the equation 



— — = cos sm (7) + - —— sm = 0. 

d(f) 2 Ha 

As before, when there is no magnetic field, H\\ = 0, the maxima occur at 
0max = ±f , ^ respectively, while the minima occur at (prain = 0, vr. 

Now, it is important to see that the search for the zeroes of the inertia 
function gives us the possibility of characterizing, at least qualitatively, the 
conditions under which tunneling can occur. Thus, starting from Eq. ()17|l . 
we verify that the zeroes of 

cos^ + -t^Hh cos = 

Ha 

will determine the situations for which tunneling will not occur since then the 
effective mass goes to infinity, and the spin will be trapped in the angular 
region around the minima of the potential. It is then direct to see that 

the roots of that equation will be = ±|,^, and = arccos ^— 

The existence of these zeroes directly imply that tunneling cannot occur; 
obviously this conclusion is verified by the exact diagonalization of the model 
Hamiltonian for the molecule, i.e., the existence of zeroes in / (0) is the 
way through which the present description reflects the diagonal behaviour of 
the phenomenological Hamiltonian in the \Sm) basis. In other words, the 
spectrum then presents a collection of doubly-degenerate energy levels. 



4.2 Fe8 cluster 

By choosing now B = —D and G = E/2, with A = 0, we have the proposed 
Hamiltonian for the Fe8 cluster ^Tl 12^1 in the absence of external mag- 
netic fields. Thus, by taking into account that, again, j = S = 10, we can 
use all the approximations discussed before so that we obtain 

Vf, (0) = -{D- E)S{S + 1) cos^ - ES{S + 1), (18) 

and 

= -TrFJ^ = -2 (^ - E) cos^ - AE, (19) 
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respectively. In the same way as in the Mnl2-acetate, we also have Vpe (0) oc 
S"^, and Ipe (0) thus also showing the relative dominance of the potential 
energy contribution to the approximate Hamiltonian. In Figure (7) we depict 
the potential and the effective mass functions for D/ks = 0.275 K and 
E/kB = 0.046 K. 

The numerical solutions of the Schrodinger equation written with those 
expressions and constants also give results that agree quite well with the 
ones obtained from the diagonalization of the corresponding original model 
Hamiltonian. The errors are less than 1% for the energy values in the lower 
part of the spectrum. 

As in the Mnl2-acetate molecule without paralel magnetic fields, the 
minima of the potential occur at = 0, vr, while the maxima occur at = 
±7r/2,37r/2. As such, the energy height defined by 

V^fT (0 = f ) - ^fT (0 = 0) = (D - E)SiS + 1) 

is obviously only dependent on the internal parameters of the molecule. 

It is now possible to evaluate the energy barrier height associated with 
the lowest energy level. If we consider the energy of the ground state (but 
not the botton of the potential well) then, 

h = -ESiS + 1) - Egs] 

using Egs = —27.645 K, obtained from the Schrodinger equation, we get hb = 
22.58 K, which is only 1.7% higher than the experimental result, 22.2 K, 
presented in fH] . 

In addition, it is important to observe that the top of the potential barrier 
lies in the negative energy domain so that there appear some energy levels - 
all those above the barrier - that seem not to be associated with tunneling 
processes. In fact there are some pairs of levels that are almost degenerate 
although this is not due to tunneling, as is also the case in the Lipkin model 
when X > Xc- This effect is due to the particular form of the Fe8 Hamilto- 
nian which embodies a subtle interplay between the term, which tends to 
produce degenerate doublets, and the (J^ — Jy) = | (J^ + J^) term which 
tends to symmetrize the energy spectrum about zero. 

In what concerns the expression that gives the inertia function, Eq.()19p. it 
is immediately seen that it presents its minima precisely at the same angles 
as the potential function does. Furthermore, for the specific values of the 
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Fe8 cluster parameters, the mass function does not exhibit infinities so that 
tunneUng is always possible, even in the absence of an external transverse 
magnetic field. In fact, the tunneling probability rate is strongly dependent 
on the behavior of the effective mass function in the barrier region 

5 Conclusions 

In this contribution we have shown that a Hamiltonian written in terms of 
an angle variable can be constructed that can describe spin systems in a 
consistent way. The construction scheme starts from a Hamiltonian given in 
terms of the operators obeying the angular momentum commutation relations 
and that is, in general, introduced from phenomenological considerations. 
The quantum Hamiltonian here proposed extends a previous description, that 
only took into account a potential energy term, by also introducing a kinetic 
energy term in which an effective mass associated with the spin system is 
present. Although this Hamiltonian has been obtained in a general form, the 
validity of the basic hypothesis supporting the assumed dominance of the first 
two terms of a moment expansion of the associated matrix energy distribution 
has been tested. Using the soluble Lipkin quasi-spin model as a testing 
ground, we have then shown that a criterion for the obtained two-terms 
Hamiltonian to be a good description of the system can be established. From 
this study it is clear that for systems with j — S > 5 the energy spectrum - 
then obtained from a Schrodinger equation - is a very good approximation to 
the one obtained from the diagonalization of the starting model Hamiltonian. 
It is then important to verify that some recently produced molecules are 
physical systems of interest that have global angular momentum values in 
the range that allow us to use the present approach. In this perspective, 
in what concerns the energy spectra and barriers, we have shown that a 
consistent description were obtained for the Mnl2-acetate as well as for the 
Fe8 cluster molecules. Being of quantum nature, this approach extends the 
semiclassical estimates and introduce essential corrections that would not be 
present. 

In what concerns the possibility of spin tunneling, we have shown that in 
the present description the potential energy function clearly exhibits energy 
barriers whose heigths are in agreement with the experimental results. At 
the same time we have emphasized the essential role played by the effective 
mass function in the discussion of tunneling in the sense that its values in 
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the barrier region are a measure of how easy it is for the system to change 
its angular orientation. The presence of infinities in the effective mass func- 
tion is an unequivocal feature indicating that tunneling is then forbidden. In 
the case of the Mnl2-acetate molecule the divergences always occur if only 
a paralel magnetic field is present, as it should be. The expression for the 
effective mass then shows the way through which the final effective Hamilto- 
nian embodies that well-known basic quantum mechanical result. The same 
is not true for the Fe8 cluster when tunneling can occur, even in the absence 
of a transverse magnetic field. The phenomenological Hamiltonian refiects 
the symmetries of the molecule, and it is precisely its particular form, with 
its coefficients obtained from the experiments, which gives rise to an effective 
mass function that is nowhere divergent. 

In the present paper we have carried out a qualitative discussion point- 
ing to the main aspects of our approach that can be used to describe the 
spin tunneling process. In a forthcoming contribution we intend to present 
quantitative estimates of the energy levels splitting associated with the spin 
tunneling. 
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Figure Captions 



Figure 1: Comparison between tlie large A^^ potential function (continu- 
ous line) and the full expression, Eq.(8), (dashed line) for the Lipkin model 
for different values of A^^. The interaction strength is x = 1-5 in all cases. 

Figure 2: Comparison between the large Ns function —I{(f)) (continuous 
line) and the full expression, Eq.(9), (dashed line) for the Lipkin model for 
different values of Ng. The interaction strength is x = 1-5 in all cases. 

Figure 3: Relative error in the ground state energies obtained from Eq. (jl3|l 
for the Lipkin model when compared with the exact results, as a func- 
tion of the interaction strength x- The continuous curve is associated with 
j = 5 (Ns = 10), the dashed curve with j = 10 (A^^ = 20), and the dot- 
dashed curve with j = 20 {Ng = 40). The lines only guide the eyes. 

Figure 4: Relative error in the ground state energies obtained from Eq. p3p 
for the Lipkin model, when compared with the exact results, as a function 
of the number of spins, Ng (= 2j), for x = 1-0. For systems with A^^ > 10 
the deviation lies below 1%. The line only guides the eyes. 

Figure 5: Curves representing the potential energy of the Mnl2-acetate 
for some different values of the paralel magnetic field. 

Figure 6: Curves representing the function I (0) of the Mnl2-acetate for 
some different values of the paralel magnetic field. 

Figure 7: Potential and effective mass functions associated with the Fe8 
cluster. 
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